Capillary Condensation and Interface Structure of a Model Colloid-Polymer Mixture 

in a Porous Medium 
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We consider the Asakura-Oosawa model of hard sphere colloids and ideal polymers in contact 
with a porous matrix modeled by immobilized configurations of hard spheres. For this ternary 
mixture a fundamental measure density functional theory is employed, where the matrix particles 
are quenched and the colloids and polymers are annealed, i.e. allowed to equilibrate. We study 
capillary condensation of the mixture in a tiny sample of matrix as well as demixing and the fluid- 
fluid interface inside a bulk matrix. Density profiles normal to the interface and surface tensions 
are calculated and compared to the case without matrix. Two kinds of matrices are considered: (i) 
colloid-sized matrix particles at low packing fractions and (ii) large matrix particles at high packing 
fractions. These two cases show fundamentally different behavior and should both be experimentally 
realizable. Furthermore, we argue that capillary condensation of a colloidal suspension could be 
experimentally accessible. We find that in case (ii), even at high packing fractions, the main effect 
of the matrix is to exclude volume and, to high accuracy, the results can be mapped onto those of 
the same system without matrix via a simple rescaling. 
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I. INTRODUCTION 



Bringing a fluid in contact with a porous medium has 
a profound influence on its characteristics and phase be- 
havior P, 0- Due to abundance of surfaces and their 
necessary proximity, surface-fluid interactions as well as 
capillarity effects play a prominent role. Moreover, the 
system may be trapped in locally stable states, and its 
behavior governed by hysteresis. Apart from the above 
fundamental questions, the study of adsorbates in porous 
media is also of great interest in applied fields ranging 
from industrial and geophysical to biomedical and phar- 
maceutical systems |2j, |3j . 

Many natural porous materials are tremendously com- 
plex on a microscopic scale: irregularly shaped pores 
build a connected void space that percolates through- 
out the sample 0,0- In contrast, to facilitate system- 
atic studies, one often relies on model pores like slit- 
like, cylindrical or spherical pores (see 0, @ and Refs. 
therein). The pore is then described conveniently in 
terms of a single parameter, its size. A different class of 
idealized system makes use of immobilized arrangements 
of fluid particles (i.e. a quenched hard sphere fluid) to 
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model a porous medium (see and Refs. therein). In 
turn, this is characterized through its density and the 
size of the spheres. However, the relevant difference to 
idealized pores is the presence of random confinement. 

The study of porous media has been focused so far 
mainly on atomic liquids. In a colloidal fluid, length 
and time scales are much larger, facilitating e.g. stud- 
ies in real space and time @. We believe that the use 
of colloidal suspensions as model systems to study the 
behavior of adsorbates in porous media can be as ben- 
eficial as their use to study many other phenomena in 
condensed matter. However, the experimental challenge 
lies in constructing three-dimensional porous media suit- 
able for colloidal suspensions. 

Colloidal porous media in 2D have been prepared by 
Cruz de Leon et. al. 0, by confining a suspension of 
large colloids between parallel glass plates. Then, these 
served as a porous matrix to a fluid of smaller particles 
which remained mobile and of which they measured the 
structure and effective potentials. To our knowledge, no 
experiment similar in spirit has been performed in three 
dimensions to date. On the other hand, Kluijtmans et. 
al. constructed 3D porous glasses of silica spheres 0, 
and silica rods i but studied the dynamics of isolated 
tracer colloids in these media. Werohski et al. studied 
transport properties in porous media of glass beads [l2T |. 
Still, such glassy arrangements of spherical colloids are 
a direct candidate for porous media suitable for colloidal 
suspensions. Sediments of large and heavy colloids as 
used in Refs 0, 0, could be brought in contact 
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with a suspension of smaller density-matched (to the sol- 
vent) colloids of which the local structure could be deter- 
mined 0,0] . However, the size ratio of the two species 
is a crucial control parameter: It has to be large enough 
(> 10) such that the small particles can penetrate the 
void space, but should still be small enough such that 
no complete separation of length scales occurs. Another 
way to realize such porous media would be to use laser 
tweezers. In a binary colloid mixture of which one of the 
species possesses the same index of refraction as the sol- 
vent (via index-matching) and the other type has a higher 
index of refraction the second species could be trapped 
while the first would still remain mobile. Using multiple 
traps at random positions in space (mimicking a fluid) 
one could then realize a model porous matrix 13] . The 
advantages of this method are the accessibility of very 
low matrix packing fractions and the full control of the 
confinement. However, the number of trapped colloids 
in such setups is typically limited to the order of 100 - 
probably too little to approach real macroscopic porous 
media, but in the right regime to be able to compare 
to computer simulations, where similar numbers are ac- 
cessible. The crucial advantage of these setups over the 
use of "natural" porous media is their model character 
arising from the use of well-defined monodisperse matrix 
spheres, while these still possess the essential features 
of random confinement and a highly interconnected void 
structure. 

One prominent phenomenon that is induced by con- 
finement is capillary condensation: A liquid inside the 
porous medium is in equilibrium with its vapor outside 
the medium. In order for a substance to phase sepa- 
rate into a dense liquid and a dilute gas phase a suffi- 
ciently long-ranged and sufficiently strong attraction be- 
tween the constituting particles is necessary. It is well- 
known that the addition of non-adsorbing polymers to 
a colloidal dispersions induces an effective attraction be- 
tween the colloids. The polymer coils are depleted from 
a shell around each colloid and overlap of these (deple- 
tion) shells generates more available volume to the poly- 
mers yielding an effective attraction between the colloids. 
Consequently, these colloid-polymer mixtures may sepa- 
rate into a colloid-poor (gas) phase and a colloid-rich 
(liquid) fluid 0. " 

The most simplistic theoretical model that has been 
applied for the study of such colloid-polymer mixture 
is the Asakura-Oosawa (AO) model 0, fl6L Il7j that 
takes the colloids to be hard spheres and the poly- 
mers to be ideal spheres that are excluded from the 
colloids. The bulk phase behavior of this model was 
studied with a variety of techniques, like effective po- 
tentials 0, [l^, free volume theory [2(j, density func- 
tional theory (DFT) [UHI and simulations [3 IM EH . 
Recent work has also been devoted to inhomogeneous 
situati ons, i.e. the free interface between demixed fluid 
phases [25|, |2(| I27L |28| , the adsorption behavior at a hard 
wall, where in particular a novel type of entropic wetting 
was found [2J, [27], |2g and the behavior in spatially peri- 



odic external potentials 29] . The surface tension between 
demixed colloid-polymer systems has been measured ex- 
perimentally and established to be much lower than for 
atomic systems [3(], EI US EH ■ Further, recent experi- 
ments confirm wetting of the colloid-rich liquid at a hard 

waii HIE!. 

DFT [36| can be used in two ways to treat adsorbates 
in porous media. The first is the (conceptually) straight- 
forward approach via treating the porous medium as an 
external potential (see e.g. Refs. [37l l38l |3!|) and to 
solve for the one-body density distributions of the fluid 
species. Those can be complicated spatial distributions, 
hence this approach is computationally demanding, but 
also yields information on out-of-equilibrium behavior 
like hysteresis in ad- and desorption curves |40L l4ll |42| . 

A recently proposed alternative is to describe the 
quenched component on the level of its one-body den- 
sity distribution |43j |. Following the fundamental mea- 
sure theory (FMT) of hard spheres IM 0, an 
explicit scheme was obtained to generate an approxi- 
mate excess free energy for (not necessarily additive) 
hard-sphere mixtures in contact with hard-sphere ma- 
trices 43]. Applied to the AO model, the results were 
compared with those from solving the so-called replica 
Ornstein-Zernike (ROZ) equations [ItL liM l49l l5(| and 
found to be in good agreement [51|. Meanwhile, this 
quenched-annealed (QA) DFT has been compared to 
computer simulations [521 ] and extended to hard-rod ma- 
trices and lattice fluids p3 |5f| . FMT in combination 
with mean field theory has also been applied to fluids in- 
side model pores |5(| |52| • 

In this article, we revisit the AO model in contact with 
a hard-sphere matrix using the QA DFT of Refs. [43Tl5l| . 
We study capillary condensation in a tiny sample of ma- 
trix as well as the fluid-fluid interface inside a bulk ma- 
trix. For both these phenomena, we distinguish two cases 
of matrices: (i) matrix particles having the same size as 
the colloids and (ii) where they are much larger. These 
correspond to the two possible experimental realizations 
we discussed earlier in the introduction, but also serve as 
representative cases because their behavior is fundamen- 
tally different. Concerning capillary condensation, we 
focus on the possible experimental realization and con- 
sider a bulk mixture in contact with in tiny sample of 
matrix. Furthermore, we elaborate if and how capillary 
condensation could be observable in such experiments. 
Concerning the fluid-fluid interface, we study the inter- 
facial profiles as well as the surface tensions inside the 
matrix. For the case of small matrix particles (i), we de- 
termine the nature of decay (monotonic or periodic) of 
the interfaces which we compare with the bulk pair cor- 
relations. For the case of large matrix particles (ii), we 
observe a simple rescaling of the bulk as well as the inter- 
face results with respect to the case without matrix. In- 
homogeneous situations like the fluid-fluid interface, are 
treated within QA DFT in a direct fashion, in contrast 
to e.g. the ROZ equations. Fluid-fluid interfaces have 
been studied before in Lennard- Jones systems in contact 
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matrix particle 



FIG. 1: Sketch of the ternary mixture of mobile colloids 
(dark), mobile polymers (transparent) and immobile matrix 
particles (grey). The polymer coils can freely overlap. There 
are three model parameters, i.e. the packing fraction of matrix 
particles, r} m and two size ratios, q = R p /R c and s = R m /Rc 
where Ri are is the radius of particles of species i. The pack- 
ing fractions of colloids and polymers, rj c and r} v respectively, 
are the thermodynamic parameters. 



with porous media using the Born- Green- Yvon equation 
as well as computer simulations (5^, and we briefly 
compare to results of our profiles. 

The paper is organized as follows. In Sec. II we define 
our theoretical model explicitly. The QA DFT approach 
is reviewed in Sec. Ill, and the results are presented in 
Sec. IV. We first consider capillary condensation in a tiny 
sample and then demixing, the interfacial profiles and 
tensions inside a matrix. We conclude with a discussion 
in Sec. V. 



II. MODEL 

We consider a three-component mixture of colloids (de- 
noted by c), polymers (p) and immobile matrix particles 
(m). Each of these particles are spherical objects with 
radii Ri and i — c,p,m and corresponding number densi- 
ties pi — Ni/V where Ni is the total number of molecules 
of species i and V is the system volume. All of these 
components are modeled as hard bodies meaning they 
cannot overlap but otherwise do not interact with each 
other, except for the polymer-polymer interaction, which 
is taken to be ideal, see Fig. ^ Consequently when r is 
the mutual distance, the pair potentials become 



Uij (r) 



oo if r < Ri + Rj 
if r > Ri + Rj 

for i, j — c,p,m except for i = j = p (1) 



and concerning the polymer-polymer interaction, this 
simply becomes 



u pp (r) = for all r 



(2) 



As all interactions are either hard-core or ideal, the 
(phase) behaviour is governed by entropic (packing) ef- 
fects and the temperature T does not play a role. The 
only thermodynamic parameters are the colloid and 



polymer packing fractions r/ c — AirR^Pc/S and r/ p = 
47ri?pP p /3, respectively. The remaining model parame- 
ters are two size ratios q — R p /R c and s = R m /R c and 
the packing fraction of matrix particles r\ m — 4irR m p m /3. 
It has to be mentioned that due to the fact that the poly- 
mers can freely overlap, the "polymer packing fraction" 
can easily be larger than one (Fig. ^| . The mixture of 
hard spheres with these last-mentioned ideal polymers 
(i.e. without the matrix particles) is called the Asakura- 
Oosawa (AO) mixture 
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III. DENSITY FUNCTIONAL THEORY 

A. Zero-Dimensional Limit 

In this subsection we derive the zero-dimensional (OD) 
Helmholtz free energy for the three-component system 
of the AO colloid-polymer mixture in contact with 
quenched hard spheres. This OD free energy is used as 
an input to construct the fundamental measure theory in 
the next subsection Here, we give only a brief derivation, 
a more extensive version with more comments is given in 
Refs. 0, EH • The essential ingredient is that we need 
to perform the so-called "double average" which refers 
to the statistical average over all fluid configurations and 
subsequently over all matrix realizations. To that end 
we consider a OD cavity which either does or does not 
contain a matrix particle. Hence, the OD partition sum 
is that of a simple hard sphere fluid, 



1 



(3) 



where z m = £ ex P(/^Mm) is the fugacity of the hard 
spheres. Further, [3 = l/k^T with fee Boltzmann's con- 
stant and fi m the chemical potential. The irrelevant pref- 
actor £ scales with the vanishing volume of the cavity but 
has no effect to the final free energy and will not be dis- 
cussed further. In general, we use an overbar to refer 
to quantities of OD systems. With the grand potential, 
(3Q m = — lnS m , the average number of matrix particles 

is T]m — Zm^P^m/ — ^m/(l ~t~ Z-m ) • 

Next, we consider the colloid-polymer mixture in con- 
tact with the matrix in zero dimensions. If the cavity is 
occupied by a matrix particle, no colloid or polymer can 
be present. On the other hand, if there is no matrix par- 
ticle, it can either be empty, occupied by a single colloid 
or an arbitrary number of polymers. Hence, 



1 matrix particle in cavity 

z c + exp(fp) no matrix particle in cavity 



(4) 



where z c and the colloid and polymer fugacities re- 

spectively. Then, the contribution — In H to the grand po- 
tential should contain the appropriate statistical weight 
for each of the cases, i.e. z m /E m for the first and 1/S m 
for the second, 



0SI 



ln(z c + exp(zp)) 
1 ~t~ z m 



(5) 
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Average particle numbers are again readily obtained via 
fji = —Zidf3£l/dzi for i — c,p (not for m). The Hclmholtz 
free energy can then be calculated using a standard 
Legendre transformation, (3F — /3f2 + J2i= c P Vi m (^i) 
and we obtain for the excess part, f3F cxc — [3F — 

PF exc (fj c , fj p ; fj m ) = (1 - fj c - fj p - fj m ) ln(l - fj c - fj m ) 
+ Vc - (1 - ^m)ln(l - fj m ). (6) 

This result can be shown to be equal from that which 
would be obtained using the so-called "replica trick" • 



B. Fundamental Measure Theory 

Fundamental measure theory (FMT) is a nonlocal den- 
sity functional theory, in which the excess part of the 
three-dimensional free energy _F CXC is expressed as a spa- 
tial integral over the free energy density 



/3F exc [{ Pl (r)}] = / dr*({nt(r)}). 



(7) 



This free energy density in turn is assumed to depend on 
the full set of weighted densities {nj,(r)}, 



<(r) = J dr'<(r-r')Mr'), 

= K®pi)0) 



(8) 



which are convolutions (denoted with <8>) with the single- 
particle distribution functions pi(r) for species i = 
c,p,m. The weight functions are obtained from the low- 
density limit where the virial series has to be recovered, 

wi(r) = e(Ri-r) 

w\ (r) = S(Ri - r) w* a (r) = 5(Ri - r)v/r 

w\{y) = 6(Ri - r)/(47rr) w^(r) = 6{R t - r)r/{Anr 2 ) 
w l (r)^S{R l -r)/(47:r 2 ) 

(9) 

with again i — c,p,m being one of the three components 
9 the Heaviside function and 8 the Dirac delta function. 
There are four scalar weight functions, with 3 to 0, corre- 
sponding respectively to the volume of the particles, the 
surface area, the mean curvature and the Euler character- 
istic and these are the so-called "fundamental measures" 
of the sphere. The two weights on the right-hand side 
of Eq. EJare vector quantities. Often a seventh tensorial 
weight is used in the context of freezing but this will not 
be used here 0, I22I |46| . The dimensions of the weight 
functions w l v are (length)' 7-3 . 

Then, the sole approximation made is that $ is taken 
to be a function of the weighted densities n l u (r) whereas 
most generally one would expect this to be a functional 
dependence. This approximation totally sets the form of 
$ and following Refs. [22LI43L l5l| we give the expression 



for $ = <f>i + $2 + &3 in terms of the zero-dimensional 
free energy derived in the previous subsection, 



$ 1= £ nMni}), (10) 

i—c,p,m 

$ 2 = Yl (n>i-<i-<2)<Pv({4}), (11) 



(12) 



with 



*>*!..*({%}) = d t 0F eKC ({fj j })/dfj i ,...dfj it . (13) 

All ^^...i-t of which more than one indices equal p are 
zero due to the form of (3F cyic . Together, Eqs. © to ifH3|) 
constitute the excess free energy functional for this QA 
system. 



C. Minimization 

Having constructed the excess free energy, we can now 
immediately move on to the the grand-canonical free en- 
ergy functional of the colloid-polymer mixture in contact 
with a matrix, 

^[Pc(r),/9p(r);p m (r)] = F cxc [p c (r), p p (r); p m (r)} + 

+ k B T J dr M r ) I ln (/^( r ) A *) - 1] + 

i—c.p 

+ E / *7>i(r)M(r)-w] (14) 

i— c,p 

Here, Aj is the "thermal volume" which is the product 
of the relevant de Broglie wavelengths of the particles of 
species i. Further, pi is the chemical potential and Vi is 
the (external) potential acting on component i. In this 
paper, we study bulk phase behaviour and the free fluid- 
fluid interfaces so we use Vi — 0. The equilibrium profiles 
are the ones that minimize the functional, 



SCI 
5pc(r) 



and 



5Q 
5p p {r) 







(15) 



This yields the Euler-Lagrange or stationarity equations 
(i = c,p), 



Pi(r) = z.exp cf ] ({pj(r)}) 



(16) 



with Zi — A i 1 exp[(3pi] the fugacity of component i and 
the one-particle direct correlation functions given by 



c, (1) (r) = 



8F oxc [{p 3 (r)}} 



= -E(S®^J (r) 

(17) 
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Obviously, the functional is not minimized with respect 
to the matrix distribution p m (r) as this serves as an input 
profile. In principle, as we are dealing with a quenched- 
annealed system in which the matrix is initially (before 
quenching) a hard sphere fluid, p ™, (r ) should still min- 
imize the hard sphere functional 0, However, as 
density functional theory allows us to generate any dis- 
tribution p m (r) by applying any suited external potential 
(which we can then remove after quenching), we do not 
need to go into the scheme of generating matrix profiles. 
Moreover, in the present paper we use fluid distributions 
of the matrix particles which minimize (at least locally) 
the hard sphere functional without external potential for 
any packing fraction. 



IV. RESULTS 

In this section, we show results of the effect of the 
hard sphere matrix on the AO colloid-polymer mixture 
concerning capillary condensation in a tiny sample of ma- 
trix, and phase behaviour and free fluid-fluid interfaces 
inside a bulk matrix. Throughout this section, we distin- 
guish between the colloid-sized matrix particles (s = 1) 
and the large matrix particles (we use s = 50). In the 
first case, as we will see, one is limited to small matrix 
packing fractions (up to r/ m of the order of 0.2) as for 
high packing fractions the pores become too small for 
the colloids and the polymers to constitute a real fluid in 
the matrix. For the large matrix particles, higher matrix 
packing fractions are accessible (up to r) m = 0.5). Fi- 
nally, in all cases we use q = 0.6, for which size ratio the 
AO model has a stable fluid-fluid demixing area (with 
respect to freezing, which we do not consider). 



A. Bulk Fluid Free Energy 

In the fluid phase, the densities are spatially homoge- 
neous, and the constant distributions Pi{r) — pt solve 
the stationarity Eqs. and i|17|) . Therefore, we only 
have to integrate the weights over space, J drwl, and the 
weighted densities become 

n\ = Vi, 

n 2 = 3r)i/Ri, 

nj = 3^/(471^), (18) 
4=3^/(4^?), 

n v2 = n vl = 0; 

with i — c,p,m. Substituting these expressions in the 
free energy density, Eq. I|10|) to (|12[1 we obtain an analyt- 
ical expression for the bulk excess free energy. Defining 
the dimensionless bulk free energy density, / = f3FV c /V 



with V c = 47ri? 3 /3 the volume of a colloid, this becomes 

77 

f(VcV P ;Vm) = Vc (hi??c - !) + -3 ( lnr /p - !) 

+ foiVcVm) ~ ^fhia(?7 c ,?7 m ) (19) 

with the last two terms being the excess free energy. 
We have separated the excess free energy in two terms 
where fo the excess free energy density of a fluid of hard 
spheres in contact with a hard sphere matrix and a is the 
fraction of free volume for the polymers in the presence 
the hard sphere colloidal fluid and the hard sphere ma- 
trix |2fj l22j |. The expressions for fo and a are quite 
extensive and given in the appendix. In going from 
Eq. to l|19l) we have discarded two terms, r\ c \n(A c /V c ) 
and (jip/q 3 ) ln(A p /V p ), linear in the colloid and polymer 
packing fractions. These have no effect on the phase be- 
haviour. Due to the ideal interactions of the polymers, 
the excess free energy density is only linear in rj p and the 
polymer fugacity becomes simply 

z pV P = yp/<x(vc,Vm)- (20) 

This relation is trivially invertible, so switching from 
system representation (using f(r/ c , rj p \ r) m )) to the poly- 
mer reservoir representation (in terms of u>(r) c , z p ;rj m ) = 
f — p, p r] p /q 3 ) is straightforwardly done. Moreover, for 
zero packing fractions of colloids and matrix particles, 
the polymer free volume fraction is trivial, a(0,0) = 1. 
Consequently, the fugacity equals the packing fraction of 
polymers in the polymer reservoir, z p V p — rj Pjr (where 
there are no colloids and matrix particles), and often, we 
use r) p ^ r when referring to the fugacity. Finally, we men- 
tion that in the absence of matrix particles, ?y m = 0, this 
theory is equivalent to the free- volume theory for the AO 
model HUHm. 

Concerning the fluid-fluid demixing, the spinodals are 
calculated in the canonical representation, by solving 
det[d 2 f / drjidrjj] = with i,j = c,p, which can be 
done analytically. Binodals are determined by construct- 
ing the common tangents of the semi-grand potential 
^(VcVpyiVm) at fixed fugacity i] p . r . 

When the matrix particles are very large, it is ex- 
pected that the excluded volume effects dominate over 
other (surface or capillary) effects. In particular, if one 
considers only one infinitely large particle, still corre- 
sponding to a nonzero matrix packing fraction, one would 
expect normal bulk behavior of the mixture as most 
of the mixture is "far" away from the matrix particle. 
Equivalently, for very large matrix particles, the total 
volume of the surrounding depletion layers, which are 
responsible for the surface effects, compared to the ac- 
tual volume occupied by the matrix particles scales with 
(47r(i? m + i? c ) 3 p m /3 — rj m )/rj m oc 3/s for the colloids and 
(47r(i? m + R p ) 3 p m /3 - Tjm) /rim oc 3q/s for the polymers, 
and these both go to zero for s — » 00. However, in this 
limit, we still need to correct for the volume as this is 
partly occupied by infinitely large matrix particles, i.e. 
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FIG. 2: Phase diagram of a bulk AO colloid-polymer mixture 
(q — 0.6) in chemical contact with a tiny sample of porous 
material (s = 1). Concerning the bulk mixture: the thick 
full curve is the binodal, the dashed is the spinodal, the large 
filled circle (where they meet) is the critical point and the 
straight (thin) lines are the tie lines connecting coexisting 
state points. The capillary lines (full curves) appear in the up- 
per left (colloid-poor) part of the phase diagram and each ter- 
minates in a capillary critical point (small filled circles). From 
lower-right to upper-left, the curves correspond to an increas- 
ing packing fraction of the matrix, rj m = 0.05,0.1,0.15,0.2 
(the last one is practically on the vertical axis near r\ v = 0.8). 
The inset shows a sketch of the setup: a test tube is sealed at 
the top and filled with the colloid-polymer mixture (densities 
r/ c and r) p ) and at the bottom lies the tiny sample of porous 
matrix. 

V — > (1 — i] m )V. Indeed, applying s — > oo to the bulk 
free energy of Eq. (|19fl . we re-obtain the bulk behaviour 
of the plain AO colloid-polymer mixture without matrix, 
i.e. it can be shown that 

lim f(r] c ,r) p ;r] m ) = (1 - T) m )f [ - — - — , — ^ — ;0 J , 

(21) 

where the free energy density has to be rescaled as well. 
This term can be considered to be the zeroth in a 1/s- 
expansion of the free energy of which higher order terms 
should correspond to effects due to surfaces, capillarity, 
curvature etc. However, because of the formidable form 
of the free energy it is a daunting task to connect every 
term to a certain phenomenon and we leave this to future 
investigation. It is worth mentioning that a power series 
in 1/s is only a simple model dependence. In general, 
there can be non-analyticities, e.g. arising from wetting 
phenomena around curved surfaces [6Cj of matrix parti- 
cles. 



B. Capillary Condensation in a Tiny Sample of 
Porous Matrix 

A porous matrix of quenched hard spheres stabilizes 
the colloid-rich phase with respect to the colloidal gas 




FIG. 3: Same as in Fig. [5] but now for a sample of porous ma- 
trix with s = 50. Here the matrix packing fractions increase 
as follows: rjm = 0.1,0.2,0.3,0.4,0.5 (again from lower-right 
to upper- left). The inset is the magnification of the area with 
the capillary critical points. 

phase [Hi] . This is called capillary condensation and it is 
due to the attractive depletion potential between the col- 
loids, which also acts between colloids and matrix parti- 
cles. In this subsection we present capillary condensation 
in a representation which is appropriate to compare with 
experiments. In an experimental setup one typically has 
a canonical ensemble, i.e. a test tube, of colloid-polymer 
mixture. By adding a tiny sample of porous material, the 
bulk mixture in the test tube acts as a colloid-polymer 
reservoir to the sample, but vice versa, if the sample is 
small enough, its state will not have any effect on that 
of the bulk mixture (sec Fig. |2finset)). In the colloid- 
poor (and polymer-rich) part of the phase diagram, on 
approaching bulk coexistence, the conditions for coexis- 
tence in the porous sample are reached before those in 
bulk, i.e. capillary condensation in the sample occurs. 
Hence, capillary condensation appears as a line in the 
system representation terminating in a capillary critical 
point. This is shown in Figs. [2] and for the case of s = 1 
and s = 50 respectively and various densities of the ma- 
trix. The coexistence of the bulk colloid-polymer mix- 
ture appears in the usual system representation, where 
tie lines connect coexisting states. For each of the matrix 
densities, a capillary line runs along the bulk binodal in 
the colloid-poor part of the phase diagram. 

First, we determine the conditions for coexistence in- 
side the matrix, i.e. we compute the combinations of 
chemical potentials /i^ ™ and A*p °oex ' f° r which demix- 
ing occurs within the porous sample. These are fixed by 
the chemical potentials of the bulk colloid-polymer mix- 
ture, /^c Ulk and /iip Ulk , so solving 

£ ulk (VcV P )=^°™, (22) 

bulk/ \ porous 

H-p V'/ci 'lp J fpjcocx 

for f] c and rj p , we obtain the capillary lines in the phase 
diagram in system representation. The trend can be 
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0.4- 
0.2 

qI , 1 , 1 , 1 , 1 , 

u 0.1 0.2 0.3 0.4 0.5 

Vc 

FIG. 4: Fluid-fluid binodals of an AO colloid-polymer mix- 
ture (q = 0.6) inside a bulk porous matrix (s = 1). Tie 
lines connecting coexisting state point are not drawn but 
run horizontal. The lower (thicker) curve is the result in 
the absence of any matrix, r) m = 0. For the other curves, 
the matrix packing fraction increases from bottom to top, 
r) m = 0,0.05,0.1,0.15,0.2. The filled circles are the critical 
points (large, r\ m — 0). The dotted line is the Fisher-Widom 
line for rj m = below of which the decay of correlations in the 
fluid is oscillatory and above these are monotonic. The point 
where the FW line hits the binodal is marked by a (large) 
star. The FW lines for the other matrix packing fractions 
are not shown, only their crossings with the binodals (small 
stars) . 



spotted from Figs. |21 and [31 increasing the matrix pack- 
ing fraction in the sample, the capillary line moves away 
from the bulk binodal but at the same time the capillary 
critical point shifts away from the bulk critical point. 
Qualitatively this applies to both the s = 1- and s = 50- 
cases. However, in the s — 50-case the capillary lines 
extend to much closer to the bulk critical point, but they 
are hardly distinguishable from the bulk binodal. Con- 
cerning the colloid-sized matrix particles, s = 1, these 
capillary lines are well separated from the bulk binodal, 
but the capillary critical points are located much deeper 
into the colloidal gas regime. 

Next, we briefly discuss the implications this has for 
possible experiments. Focusing on the case of the large 
matrix particles, s = 50, we take as an example r\ m = 0.5. 
In this case the difference in chemical potential at co- 
existence of the mixture in bulk and inside the porous 
sample at constant polymer fugacity is of the order, 
/3A/i™ ex ~ 0.1, and it scales roughly with l/s. This 
difference is very small and brings up the question if 
this (i.e. capillary condensation) is observable in exper- 
iments. Typically, the effect of gravity is reduced by 
density-matching the colloids with the solvent, i.e. can- 
celing gravity by means of buoyancy. However, this 
density- matching is never perfect, and the length scale 
(/3m c <7) _1 is a measure for its success (at infinity it is 
perfect). Here g is the gravitational acceleration and 
m c = (pcoiloid - Psolvent)^c the effective mass of the colloid 




4) 0.1 0.2 0.3 0.4 0.5 
A I , 1 , 1 , 1 , 1 , 1 

u 0.1 0.2 0.3 0.4 0.5 

Vc 

FIG. 5: Same as in Fig. |U but now for s = 50. Ma- 
trix packing fractions increase from right to left, r\ m = 
0,0.1,0.2,0.3,0.4,0.5. Inset: same curves now rescaled, i.e. 

T] p , r VS. 7J C /(1 - rjm)- 



in solution, with p C oiioid and p S oivent the mass densities 
inside the colloid and of the solvent respectively. There- 
fore this length scale is strongly dependent on the colloid 
size, (j3m c g) oc R~ 3 , and can range from microme- 
ters (large colloids) to meters (small colloids) in experi- 
ments |32j . Typically, polymers are much less sensitive to 
gravity as long as the solvent is good. When there is co- 
existence inside a test tube, there is only real coexistence 
at the liquid-gas interface whereas below and above the 
colloids have slightly different chemical potentials due to 
their gravitational energy. Consequently, moving upward 
from the interface, say Az, the colloid chemical potential 
is ([3m c g)Az lower than at coexistence. By placing the 
porous sample within Az* = /3A/i™ ex / (f3m c g) of the in- 
terface, capillary condensation should take place. Taking 
as an example, /3A/x™ ex ~ 0.1 and (/3m c g) _1 ~ 1 meter, 
it becomes clear that, within the context of this (ideal- 
ized) model, values of Az* ~ 0.1 meter should be accessi- 
ble in experiments, meaning that capillary condensation 
could in principle be observed. Complete wetting of the 
large matrix spheres by the colloidal liquid may preempt 
capillary condensation close to the bulk critical point. 
Moving sufficiently far away from the critical point to- 
wards the dilute gas regime, the effects due to complete 
wetting should disappear while capillary condensation is 
retained. 



C. Phase Behaviour inside a Bulk Porous Matrix 

We now return to the full ternary mixture in bulk, i.e. 
where in the previous subsection, the matrix was only 
a tiny sample immersed in a large system of colloid- 
polymer mixture, in this and the following subsections 
we consider the colloid-polymer mixture in a system- wide 
matrix. In this subsection, we revisit the demixing phase 
behavior which we need in the next subsections where we 
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FIG. 6: Colloid density profiles (V c = ^irKf) normal to the 
free fluid-fluid interface for increasing matrix packing frac- 
tions at fixed polymer fugacity. Parameters are: q = 0.6, 
s = 1 and r] Pt r = 1 (see Fig. 2J. The matrix packing frac- 
tion increases from top (thick profile, r\ m = 0) to bottom: 
rj m = 0, 0.05, 0.1, 0.15, 0.2. Inset: corresponding polymer pro- 
files (p P Vp vs. z/Rc, with V p = ^TvRp) for the same values 
of the matrix packing fractions (also increasing from top to 
bottom). 



study the fluid-fluid interface inside a matrix. Fig. is 
the phase diagram in the polymer-reservoir representa- 
tion for colloid-sized matrix particles, s = 1, for various 
matrix densities. Increasing the matrix packing fraction, 
there is less volume available to the colloids and the crit- 
ical point shifts to smaller colloid packing fractions. At 
the same time, the porous matrix acts to keep the mix- 
ture "mixed" and therefore, the critical point shifts to 
higher polymer fugacities. For the case of s = 1, we can 
not go to much higher packing fractions than rj m ^0.2 
as then the critical fugacity shoots up dramatically to 
unphysically large values. This may be partly due to the 
relatively large depletion shells around the matrix parti- 
cles which cause the pore sizes to become too small for 
the colloids and polymers to enter the matrix. In case 
of large matrix particles (s = 50, see Fig. the lat- 
ter effect is negligible and the pore sizes are always large 
enough. Consequently, only the excluded volume remains 
and rescaling the binodals with (1 — rj m ) is very effec- 
tive practically mapping the binodals onto each other, 
Fig. [5f inset). This rescaling is unsuccessful for s = 1 as 
can be directly seen from the fact that the critical fugaci- 
ties in Fig. 0| are different for each of the matrix densities. 

In addition, we have determined the nature of the 
asymptotic decay of pair correlations of the fluid inside 
the matrix [r3l| . These can either be monotonic or peri- 
odic and the corresponding regions in the phase diagram 
are separated by the Fisher- Widom (FW) line, at which 
both types of decay are equally long-range. This line can 
be determined by studying the pole structure of the to- 
tal correlation functions h%j in Fourier space |olj . In the 
present case of QA systems, rather than using the usual 



FIG. 7: Same as in Fig. M but now for s = 50. Other param- 
eters are: q = 0.6 and rj PtT = 1 (see Fig. @). The matrix 
packing fractions increase from top (thick profile, rj m = 0) 
to bottom: rj m — 0, 0.1, 0.2, 0.3, 0.4, 0.5. Inset: magnification 
of the rescaled profiles for the same curves, i.e. p c V c / (1 — rjm) 
vs. z/Rc (where again, matrix packing fractions increase from 
top to bottom). 



Ornstein-Zernike equations, one has to use the replica- 
Ornstein-Zernike (ROZ) equations Neglecting cor- 

relations between the replicas, these are 



hij(r) = oy(r) 



E 

t—c,p,m 



i) /i mm )(r) 
pt(cit h tj )(r) 



(23) 



with i,j = c,p, m except i = j = m. Here, for ij ^ mm, 
the Cjj(r) = —5 2 F exc /Spi(r)5pj(r) are the direct correla- 
tion functions for which we obtain analytic expressions 
by differentiating Eq. [7] The matrix structure is deter- 
mined before the quench, so c mm and h mm are those 
of the normal hard sphere fluid at density p m (Percus- 
Yevick-compressibility closure, see Refs. [Ij, This 
analysis follows closely that of Ref. in which more de- 
tails are given. In view of our subsequent interface study, 
we focus on the point where the FW line meets the bin- 
odal. In Fig.0] (s = 1), these are denoted by stars, and 
we observe that the shifts due to the matrices follow the 
same trend as the critical points. In case of s — 50, we 
have not determined the FW lines, but there is no reason 
to expect the simple rescaling of the case without matrix 
to fail in this case. Furthermore, concerning the density 
profiles (in the next subsection, s = 50), we stay well 
within the oscillatory regime. 



D. Fluid-Fluid Profiles inside a Bulk Porous 
Matrix 

We have calculated density profiles at coexistence nor- 
mal to the colloidal gas-liquid interface. In this case 
of planar interfaces, the density distribution is only a 
function of one spatial coordinate z; i.e. pi(v) = Pi(z). 
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FIG. 8: Fluid-fluid surface tensions vs. the difference in colloid 
packing fractions of the two fluid phases for q — 0.6, s = 1 
and various values of the matrix packing fraction. The matrix 
packing fractions increase from right (r\ m = 0, thick curve) to 
left, r)ra — 0, 0.05, 0.1, 0.15, 0.2. Each curve is computed from 
the critical point, rj Pt ,. — rjp"^ (where 7 = 0) until twice the 
critical fugacity, rj Ptr — 2rj p " r lt \ 



FIG. 9: Same as in Fig. IHJbut now for s = 50. Again, matrix 
packing fractions increase from right (r\ m = 0, thick curve) 
to left, r\m = 0,0.1,0.2,0.3,0.4,0.5. Each curve is computed 
from the critical point, r] p , r — ?7pjr (where 7 = 0) until twice 
the critical fugacity, r/ p , r — 2r] p m r lt \ Inset: the same (but 



rescaled) curves are shown, i.e. 



pyR£/(i - nm) vs. \4 



(II) 



The only dependence on the other two degrees of free- 
dom is in the weights and this can be integrated out, 
to obtain projected weights, w l v (z) — J dxdyw\,{r) (see 
e.g. The profiles are discretized and calculated 

via an iteration procedure, i.e. we insert profiles on the 
right hand side of Eq. (|16fl and then obtain new pro- 
files on the left hand side, which are then reinserted 
on the right. Using step functions as iteration seeds, 
this procedure converges in the (local) direction of the 
lowest free energy. We normalize the densities as in 
bulk, i.e we plot pi{z)Vi so that pi{±oo)Vi — rjf' U \ 
with I and II referring to the coexisting phases. The 
zero of z is set at the location of the interface, de- 
fined through the Gibbs dividing surface of the colloids: 

Loo dz [Pc{z) - Pc(-Oo)} + f °° dz[ Pc (z) - p c (oo)] = 0. 

In Figs. El and [7| we have plotted the colloid profiles 
normal to the interface for s — 1 and s = 50, respectively. 
Colloid profiles are shown for increasing densities of the 
matrix at fixed fugacity, ry PjI . = 1, corresponding to the 
bulk binodals in Figs. 0] and [SJ For the case of s = 1, this 
means that, as the critical point shifts to higher fugaci- 
ties, the profiles are effectively taken at fugacities closer 
to the critical value. We observe this well-known behav- 
ior in Fig. HjJ close to the critical point the profiles are 
smoother and modulations less pronounced. Away from 
the critical point, the interface is sharp but the periodic 
modulations due to the surface extend to far in the bulk 
fluid. The inset of Fig. H3 shows the corresponding poly- 
mer profiles. In Ref. |58j . the main result is that the 
interface widens due to the porous medium. The same 
happens here and is due to fact that one is effectively 
closer to the critical point. 

In Fig. [71 as we saw for the bulk phase diagram, there 
is a simple rescaling at work and the profiles merely dif- 



fer with a factor (1 — 7y m ). The inset in Fig. [7| shows 
the same colloid profiles but now rescaled and we have 
zoomed in on the region close to the interface. Clearly, 
even the modulations follow the case without matrix with 
the same accuracy as the bulk coexistence values in the 
inset of Fig. [SJ 

We have also studied the asymptotic decay of correla- 
tions with the interface via the density profiles. These 
must be of the same nature as the decay of the direct 
correlations in bulk (determined via the ROZ equations, 
see previous subsection), i.e. either monotonic or peri- 
odic |6l| . However, determining the crossing points of 
the FW line with the binodals using the interfacial pro- 
files yields a systematic shift away from the critical point, 
compared to the bulk calculation (~ 5%). Probably, this 
is due to numerical limits. Close to this crossing point 
both (the periodic and the monotonic) modes of decay are 
equally strong, so only far away from the interface truly 
asymptotic behavior may be observed. However, there, 
the periodic modulations may have become too small to 
be observable. Furthermore, our numerical routine has 
no real incentive to minimize the tails of the profiles as 
the gain in free energy is very low. 



E. Fluid-Fluid Surface Tension inside a Bulk 
Porous Matrix 

The presence of the matrix also affects the surface ten- 
sion between the colloidal liquid and gas phases. The 
interfacial or surface tension 7 of planar interfaces in the 
grand canonical ensemble is defined through 

7 A = fi inh + PV (24) 
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where A is the amount of surface area, Jlinh is the grand 
potential for the inhomogeneous system and P the pres- 
sure (i.e. — PV is the grand potential for the homoge- 
neous bulk system). With our numerical scheme we cal- 
culate density profiles in z-direction so it makes sense to 
write the surface tension as an integral, 

7 = J dz [uj(z) + P] (25) 

with 

u(z) = k B T P*( z ) [MM^Ai) - 1] 

i—c,p 

- £ + fcBT$(K(*)}). (26) 

i— c,p 

The quantity lo(z) is a "local" grand potential density 
whose average over space yields the actual grand poten- 
tial per unit of volume flwh/V. In Figs. 0D and ED we 
have plotted the surface tension versus the colloidal den- 
sity difference in the two phases for s = 1 and s = 50, 
respectively. In both cases the effect of the matrix is 
that the surface tensions increase faster with the differ- 
ence l^c 11 "* — T)c\ which is of course due to the fact that 
the coexistence area becomes less wide as the coexist- 
ing packing fractions themselves become smaller. In the 
inset of Fig. [5J we show the same curves rescaled with 
(1 — Vm), and the rescaled graphs fall almost on top of 
the original one without any matrix. Here, we note that 
also the surface tension has been rescaled with (1 — T) m ); 
this is needed from Eq. (|21[1 as the free energy density 
uj(z) needs to be rescaled as well. Again, this rescaling 
procedure is not successful for s = 1. 

Often, the surface tension is plotted against the rela- 
tive distance to the critical point, (^ P ,r/%r — 1) [13 • 
However, this does not improve the rescaling for s = 1 
and this can be seen from the fact that the end points 
of the curves in Figs. Eland El are all at twice the critical 
fugacity, r/ p ^ r = 2i}p^ and the surface tensions (rescaled 
or not) are at quite different values at the end points. 

V. CONCLUSION 

We have considered the full ternary system of hard 
spheres and ideal polymers (represented by the AO 
model) in contact with a quenched hard sphere fluid act- 
ing as a porous matrix. Using a QA DFT in the spirit of 
Rosenfeld's fundamental measure approach, we studied 
capillary condensation in a tiny sample of matrix as well 
as the the fluid-fluid interface inside a bulk matrix. The 
results have been presented in terms of two types of ma- 
trices: (i) colloid-sized matrix particles (size ratio s = 1) 
being a reference system and (ii) matrix particles which 
are much larger than the colloids (size ratio s = 50). 
The case of small matrix particles is limited to relatively 
low packing fractions (rj m ~ 0.2), whereas in the second 



case, much higher matrix packing fractions are accessible 
(jlm ~ 0.5), the pores of the matrix being much larger. 
Additionally, we have suggested that case (i) as well as 
(ii) could in principle be realized experimentally in 3D, 
i.e. using laser tweezers and colloidal sediments respec- 
tively, to serve as a model porous medium for colloidal 
suspensions. 

We have shown that in the limit of infinitely large ma- 
trix particles, the standard AO results (without matrix) 
are recovered via a simple rescaling. In case of s — 50 our 
bulk but also the interface results can be mapped onto 
the case without matrix with high accuracy. However, in 
the case of small matrix particles (s = 1) this mapping 
fails, which is due to the more complex (and smaller) 
pore geometry on the colloidal scale. 

Assuming a more "experimental" point of view, we 
have considered a tiny sample of porous matrix immersed 
in a large system of colloid-polymer mixture. When the 
fluid-fluid binodal is approached in the colloid-poor re- 
gion of the phase diagram, capillary condensation occurs 
in the sample. This transition appears as a capillary line 
in the phase diagram (in system representation) extend- 
ing along the binodal and ending in a capillary critical 
point. In case of small matrix particles, the capillary 
lines (for various densities of the matrix) are well sep- 
arated from the bulk binodal but the capillary critical 
points lie deep into the colloidal gas regime. Concerning 
the large matrix particles, these capillary critical points 
are located closer to the bulk critical point, however, the 
capillary lines are also very close to the binodal. Still, us- 
ing density-matched colloidal suspensions, we argue that 
capillary condensation may be observable in experiments. 

We have computed fluid-fluid profiles inside the porous 
matrix as well as the corresponding surface tensions. For 
s = 50, these can be mapped onto the case without 
matrix but for s = 1 the critical point shifts to higher 
polymer fugacities. Therefore, increasing the density of 
the matrix, profiles become smoother due to effective ap- 
proach of the critical point. Solving the ROZ equations, 
we have also determined the crossover between mono- 
tonic and periodic decay of pair correlations of the mix- 
ture inside the matrix for s = 1. Comparing these with 
the decay of the interfacial correlations we find a small 
discrepancy which is probably due to numerical limits. 

It should be noted that we do not expect our current 
approach to satisfactory describe the (subtle) phenomena 
associated with wetting of the curved surfaces of the ma- 
trix particles by the colloidal liquid Ell ■ Especially 
for s — 50 close to the critical point in the complete 
wetting regime (of the planar hard wall), we can well 
imagine that the growths of thick wetting films preempt 
our capillary condensation transitions as well as disturb 
the fluid-fluid interfaces. 

Concerning the fluid profiles, we have only considered 
a homogeneous background of matrix particles in this pa- 
per. It would be interesting to use inhomogeneous matrix 
realizations, as e.g. a step function of zero and nonzero 
matrix packing fraction (i.e. the interface of empty space- 
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matrix) or a constant matrix background in contact with 
a hard wall. Both types could give rise to interesting 
and substantially modified wetting behaviour. Addition- 
ally, one could also consider other types of matrices, e.g. 
quenched polymers or combinations of quenched colloids 
with quenched polymers 0, E3 • These are maybe less 
realistic from an experimental point of view but still in- 
teresting due to the competition of capillary condensa- 
tion with evaporation. 

As we have mentioned in the introduction, there 
are no experiments concerning phase behavior of col- 
loidal suspensions in contact with 3D porous media to 
our knowled ge. We hope that the accumulating re- 
sults 0, H HM El d El H3, including those in this 
paper, may encourage more experimental efforts in that 
direction. It is important to keep in mind that a suitable 
porous matrix is a compromise between length scales: 
large enough to allow penetration of the colloids into the 
void space, but small enough to retain significant surface 
and capillary effects. In colloidal fluids in general, these 
last-mentioned effects are known to be much smaller than 
in atomic systems, thus providing a formidable challenge 
to experimentalists aiming to observe e.g. capillary con- 
densation of a colloidal suspension in a porous matrix. 
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APPENDIX A: BULK FLUID FREE ENERGY 

The bulk free energy of the colloid-polymer mixture in 
contact with a homogeneous hard-sphere porous matrix 
as given in Eq. (|T9"|) is 



firtdVpiVm) = Vc 0?7c - 1) + -§ (hir/p - 1) 

+ MVcVm) ~ ^lna(i] Cl r] m ) (Al) 



We note the occurence of only third and lower powers of 
1/s in both fo and a which are given by 



fo{Vc, rim) = ^ hi(l - r/m) - + ?/ c ) ln(l - r) c - ?? m ) + 



3r) c ri m (2 + ri c (r) m - 2) - 2r/ m ) 



3?7c 



2s 3 (1 - ?7 C - r/ m y 



2s 3 (l-r)m) 2 (l-Vc-Vm) 2 
{r] m (2 - 2r/ c + r]m)s + rj m (2 + rjc - 2r] m )s 2 + tj c (2 - rj c - 2r) m )s 3 } (A2) 



and 



lna(77 c ,77 m ) = ln(l - r/ c - i] m )+ 



2s 3 (l - rj c - r? m ) 3 



{277 m (l + rjl + n m + T7, 2 n - 77 c (2 + 77 m ))<7 2 



-377 m9 [-2 + rim + 7/ 2 n + 27/ c 2 (-l + q) - 77 c (-4 + r) m + 2q + Ar) m q)] 



,21 s 2 



+677 m [(1 - ?7 C - 77 m ) 2 + 3t7 c (1 - 7/c - r)m)q + i? c (l + 2?7 C - 77 m )g 2 
+77 C [-77 C (1 - 77 m )(12 + (3 - 2q)q) + 2(1 - 77 m ) 2 (3 + g(3 + q)) + r? 2 (6 + q(-3 + 2q))] s 3 } . (A3) 
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